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Abstract 

A novel expansion — which generalizes Magnus expansion — of the 
evolution operator associated with a (in general, time-dependent) per- 
turbed Hamiltonian is introduced. It is shown that it has a wide range 
of possible solutions that can be fitted according to computational con- 
venience. The time-independent and the adiabatic case are studied in 
detail. 
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1 Introduction 



The explicit determination of the evolution operator associated with a quantum 
system is a 'touchy business'. If the Hamiltonian of the system does not depend 
on time and has the form of the sum of a solvable unperturbed Hamiltonian 
plus an analytic perturbation, one can use the tools of standard perturbation 
theory PP|2] for linear operators, based on the expansion of the resolvent, in 
order to achieve approximate expressions of the cvolutor, or apply a suitable 
unitary operator perturbative approach |] ^ 

If the Hamiltonian is time-dependent (i.e. it describes a non-isolated quantum 
system), the problem is in general even more radical. In fact, it is well known 
that, whenever the values of the Hamiltonian at different times do not commute, 
the evolutor does not admit a simple formal expression. 

In two fundamental papers 6 in the history of quantum electrodynamics, 
Dyson developed an expansion of the evolution operator that has been adopted 
extensively in any field of physics. Dyson expansion has a transparent physical 
interpretation in terms of time ordered elementary processes which makes its 
application particularly appealing, especially in quantum field theory. On the 
other hand, for many applications, Dyson expansion has severe drawbacks, as a 
low convergence rate and the lack of unitarity of its truncations [7] . 
Later, Magnus [HI introduced an expansion of the evolution operator such that 
each of its truncations retains the property of being unitary. Magnus ex- 
pansion has been 'rediscovered' and re-elaborated several times (see for in- 
stance ref. [7]), applied successfully to several problems — nuclear magnetic 
resonance % , atomic collision theory |10j . molecular systems in intense laser 
fields neutrino oscillations in matter to quote just a small sample — 
and its convergence properties have been studied J3j [T3] . 

On our opinion, Magnus expansion, rather than Dyson expansion, should 
be regarded as the most natural generalization of the expression of the evolutor 
associated with a time-independent Hamiltonian. In fact, it is written in the 
form of the exponential of the expansion of a suitable time-dependent anti- 
hcrmitian operator which can be deduced, order by order, from the Hamiltonian 
of the system. Now, precisely for this reason, just like for the evolutor generated 
by a time-independent Hamiltonian, the problem of computing explicitly the 
action of (any truncation of) the Magnus expansion on the state vectors is non- 
trivial. Expanding the exponential would lead to non-unitary truncations, thus 
to the loss of the most important feature of Magnus expansion. Then, the issue 
of finding a generalization of Magnus expansion retaining the property of having 
unitary truncations, but allowing simpler explicit solutions, arises in a natural 
way. 

In the present paper, we have tried to achieve this result. Our basic idea is 
simple: to combine the Magnus expansion with the passage to a suitable interac- 
tion picture. Precisely, given a perturbed Hamiltonian, after the usual passage 
to the interaction picture that decouples the unperturbed dynamics (which is 
supposed to be known explicitly), one switches to a further interaction picture, 
depending on the perturbative parameter, in order to achieve computational 
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advantages. We stress that the idea of 'adapting' Magnus expansion is not 
new. It appears in a paper by Casas et al. |15| in which the authors introduce 
the Floquet- Magnus expansion for the evolution operator associated with a (in- 
teraction picture) Hamiltonian depending periodically on time. Our approach 
generalizes the one proposed by Casas et al. even in the case when the interac- 
tion picture Hamiltonian is periodic on time. 

We have made the choice of skipping mathematical complications. For instance, 
it is known that even a simple passage to an interaction picture can be mathe- 
matically tricky (see, for example, ref. ^H]). Our choice is motivated by various 
reasons. First of all, we believe that heuristic investigation should always pre- 
cede rigorous re-elaboration. Once that it is clear what the basic 'rules of the 
game' are, one can adopt the most appropriate mathematical tools. Moreover, 
we avoid the risk of hiding in a cloud of technicalities the main ideas and of 
discouraging those physicists who may want to apply our method for solving 
problems. It should be also observed that a recent trend in quantum mechan- 
ics is to focus on systems which can be described by effective Hamiltonians in 
finite-dimensional Hilbert spaces (consider, in particular, the huge research area 
related to quantum computation and quantum information theory; see ref. |17j 
and the rich bibliography therein) . The study of these systems is not affected by 
all the technicalities associated with the infinite-dimensional spaces but retains 
all the most intriguing features of quantum physics. 

The paper is organized as follows. In section [3 we introduce the basic 
decomposition of the evolution operator which will allow us to obtain a pertur- 
bative expansion. Two important cases — the time-independent case and the 
adiabatic case — will be considered in sections OH and In section 03 we will 
study the general time-dependent case and show how the solutions obtained 
correspond, in the time-independent case, to the ones obtained in section 

2 Basic assumptions and strategy 

Let us consider a time-dependent perturbed Hamiltonian H(X;t), namely a 
selfadjoint linear operator of the form 

H(X;t)=H (t)+H o (X;t), (1) 

where Ho(t) is a selfadjoint (and, in general, time-dependent) operator — the 
'unperturbed component' — and H&(X;t) is a time-dependent perturbation; 
precisely, we will assume that A i— » H (X]t) is (for the perturbative parameter 
A in a certain neighborhood of zero and for any t) a real analytic, selfajoint, 
bounded operator-valued function, with H o (0;t) = 0. A real analytic function 
can be extended to a domain in the complex plane. Keeping this fact in mind, 
we will specify that a given property holds for X real. For instance, the analytic 
function A i— ► H (X;i) will take values in the selfadjoint operators for A real 
only. 

Let U (A; t, to) be the evolution operator associated with H (A; t), with initial 
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time to] namely (K = 1): 

iU(X;t,t )=H(X;t)U(X;t,t ), U(X; t, t ) = Id, (2) 

where the dot denotes the time derivative. Then, we have that 

U{X;t,t ) = U {t,t )T(X;t,t ), (3) 

where Uo{t, to) and T(A; t, to) are respectively the evolution operator associated 
with the unperturbed component Ho(t) (evolution operator which, if the un- 
perturbed Hamiltonian is time- independent, Ho(t) = Ho, is obviously given by 
e -iH (t-t )^ anc j evolution operator associated with the interaction picture 
Hamiltonian 

H(X; i, t Q ) := U (t , t) H (X; t) U (t, t ). (4) 
Let us notice explicitly that, since H(0;t,to) = 0, we have: 

T(0;Mo) =W. (5) 

We will suppose that the unperturbed evolution Uo(t, to) is explicitly known. 
Then the problem is to determine perturbative expressions of T(A; t, t a ). To this 
aim, the central point of the paper is the assumption that T(X;t,to) has the 
following general form: 

T(A; t, t ) = exp (-i Z(X; t, t )) exp C(A; t, t ) dij exp (i Z(X; t , t Q )) , (6) 

where (A; t) Z(X;t, to), (X;t) t— > C(X;t,to) are operator-valued functions 
which depend analytically on the perturbative parameter A; in agreement with 
condition JSJ), we set: 

Z(0;t,t ) = 0, C(0;t,io) = 0, Vt. (7) 

We stress that the presence of the term exp(i Z(X; to, to)) in formula I© ensures 
that T(X;to,to) = Id, allowing the possibility that Z(X;to,to) ^ 0. 

It will be seen that decomposition has a wide range of solutions and that 
a possible choice for fixing a certain class of solutions is given by imposing the 
condition C(X;t,to) = C(A), i.e. assuming that the function (A; t) i— > C(X;t, to) 
does not depend on time. This decomposition includes, as particular cases, 
two decompositions of the evolution operator that have been considered in the 
literature: 

• the decomposition that is obtained setting 

Z(X;t,to) =0, Vt, 

in formula © , decomposition which is at the root of the Magnus expansion 
of the evolution operator [S]; 
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• the classical Floquet decomposition that holds in the case where the in- 
teraction picture Hamiltonian depends periodically on time (let us denote 
the period by T) — decomposition which is obtained setting 

C(X;t,t Q ) = C(X), Z(A;to,*o) = 0, 

and assuming that (A, t) t— > Z{X\ t, to) is periodic with respect to time with 
period T — and that is at the root of the Floquet-Magnus expansion of 
the evolution operator |15| . 

From this point onwards, for notational convenience, we will fix to = 0. 
Then, decomposition © can be rewritten as 

T(A; t) = exp (-iZ(X; *)) exp (~iJ Q C ( A ; l ) dt ) ex P ( iZ W) > (8) 
where we have set: 

T{X;t)=T(X;t,0), Z(X; t) = Z(X; t, 0), Z(X) = Z(X;0), C(A; t) = C(A; t, 0). 

We are now ready to obtain a perturbative expansion of T(X;t). In fact, if 
we require the interaction picture evolution operator to satisfy the Schrodinger 
equation, we get: 

H{X;t)T(X;t) = if(X;t) 

= e -iZ(X;t)J^ e isZ(X;t) z(X;t)e- lsZ ^)ds e -^C(X;t)di e iz(X) + 

+ e -^(^*^ 1 (e-^o c ( A;t ) dt C(A;t)e^o c ( A;t ) dt )d s e -^C{X;t) di &i z (X) ^ (Q) 
where we have used the remarkable formula (see, for instance, ref. |18j 1 

^ e F = e F ^ (e- sF F e sF ) ds = J* (e sF F e- sF ) ds e F ', (10) 



which extends to an operator-valued function t i— > F(t) the standard formula for 

the derivative of the exponential of an ordinary function. Next, let us apply to 

each member of equation © the operator e lZ ^ x '^ on the left and the operator 
e -iZ{X) e iC{X)t Qn tne rig]lt . we find . 



Ad oxp ( iZ ( A . t )) H(X; t) — I ^Ad cxp (j sZ ( A;t)) Z(X;t)jds 

( A d cxp (-^c(A; t )d t )^(A;i))d S , (11) 

where we recall that, given linear operators X, Y , with X invertiblc, 

Adjt-F := xyx" 1 . (12) 
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Then, since X is of the form e x , we can use the well known relation 

Ad cxp(x) y = exp(ad x )r = ^-ad^y, (13) 

fc=0 

with adx denoting the k-th power (ad^- = Id) of the superoperator adx defined 

by 

ad x Y :=[X,Y\. (14) 

Eventually, applying formula (|13|l to equation ljll|l and performing the integra- 
tions, we obtain: 

00 -k 00 -k 

E ~L1 ad l(A; t ) H(X; t) = E t^W ad* (A . t) Z(X; t) 

+ E(^T)! ad W)d t ^;^)- (15) 

This equation will be the starting point for the determination of the operator- 
valued functions (A, t) i— > Z(X; t) and (A, t) i— » C(A; i) at each perturbative order 
in A, task that will be pursued sistematically in the next sections. 

3 The time-independent case and the perturba- 
tive adiabatic approximation 

In this section, we will consider two important cases: 

1. the case where the Hamiltonian 0) does not depend on time; 

2. the case where the Hamiltonian is slowly varying with respect to time 
(so that the adiabatic approximation can be applied). 

As we will show later on, the second case, within the adiabatic approximation, 
can be treated by a method analogous to the one adopted in the first case. In 
both cases it will be convenient to set 

Z(X;t):=U (t)Z(X;t)U (t)\ Z(X) = Z(X; 0) = Z(X), (16) 

and re-express equation (|15JI in terms of the transformed operator Z(X;t). To 
this aim, let us first notice that 

Z(X; t) = kd Ua(ty (Z(A; t) - i ad 2(A;t) H ) . (17) 

Besides, given linear operators X, X and Y, with X invertible, one can show 
inductively that 

ad k AdxX Ad x Y = Ad^ad^F, k = 0,1,2,... . (18) 
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Then, since Z(X;t) — Ad Uo uy Z(X;t), H(X;t) = Ad Uo up H (X;t) and rela- 
tion 117|l holds, using formula 1181) . from equation (|15() we obtain: 

00 -k 00 -k 

M u (ty E T\ ad 2(A;*) H o( X > *) = Ad ^o(*)t E JkTTv ad ^^*) i(A; *) 

°° jfe 

- Ad C/o(t)tE^? ad 2(A;t)^0 
fc=l 

+ E(^T)!^C(A;i)dt^^)- (19) 

Next, applying the superoperator Ad[/ (t) to each member of this equation and 
rearranging the terms, we get 

00 

E \ x ad* (A;t) (tf (i) + H (X: t)) + H (X; t) = Ad Uo(t) (c(A; t) 

fc=l 

OO / . 

+ E^^(W C M 
+ E7^ a 4(A; t ) i (^). (20) 

fc=0 v 

3.1 The time-independent case 

In the time-independent case, we have that H (t) = H , ff (A;t) = i7 (A), and 
it is natural to set: 

C(X;t) = C(X), Z(X;t)=Z(X;Q) = Z(X;0) = Z(X). (21) 

Then, equation (|2U[) assumes a much simpler form: 

00 jk 

E adS (A) (ffo + H*(A)) + ff«(A) - e^ ' C(A) ^ (22) 
fc=i 

Now, observe that the first member of this equation does not depend on time, 
hence the function t 1— > er %Hot C(A) e lfl °* must be constant. It follows that, if 
we want equation 122[l to be consistent, we have to assume that 

[C(X),H ]=Q, (23) 

i.e. that C(A) is a constant of the motion for the unperturbed evolution generated 
by Hq. Eventually, we obtain the following fundamental formula: 

00 .u 

E ^ ad| (A) (ff + H (X)) + H (X) = C(X). (24) 

k=l 



7 



At this point, we are ready to obtain perturbative expansions of the operators 
C(A) and Z(X) (hence, of the interaction picture evolution operator T(A; <)). We 
will assume that the unperturbed Hamiltonian Hq has a pure point spectrum, 
while the case where this hypothesis is not satisfied is a particular case of the 
general treatment developed in section [S] (indeed, in this case the formulae 
obtained in this section, which involve eigenvalues and eigenprojectors, make no 
sense). We will denote by E\, E 2 , ■ ■ ■ the (possibly degenerate) eigenvalues of H 
and by Pi, P2, ■ ■ ■ the associated eigenprojectors. Since the functions A 1— ► P (A), 
A i-> C(A) and A >-> Z{\) are analytic and H (O) = (7(0) = Z(0) = 0, we can 
write: 

00 00 00 

P (A) = ^A"P„, C(A) = Z(\) = J2* n Zn- (25) 

n— 1 n— 1 n— 1 

Now, in order to determine the operators {C n } n gN and {Z n } nt =fi, let us sub- 
stitute the power expansions l|25|) in equation Q24[l: in correspondence to the 
various orders in the perurbative parameter A, we get the following set of con- 
ditions: 

C 1 -i[Z 1 ,H ]-H 1 =0 , [d,H Q }=0 (26) 
C 2 -i[Z 2 ,Ho] + ^[Z u [Z 1 ,H }]-i[Z 1 ,H 1 }-H 2 =0 , [C 2 ,P ] = (27) 

where we have taken into account the additional constraint [C(X),Hq] = 0. 
This infinite set of equations can be solved recursively and the solution — as it 
should be expected (we will clarify this point soon) — is not unique. The first 
equation, together with the first constraint, determines Z\ up to an operator 
commuting with Hq and C\ uniquely. Indeed, since 

[C u H a ]=0 => Cx = J2P m CxP m (28) 

m 

and 

[Z U H ]=J2 ( E l - E i ) P i Z ^ p i> ( 29 ) 

we conclude that 

ci=J2 p ™ H i p ™ ( 3 °) 

m 

and 

Zi=J2 Pm Zl Pm + { J2( El - P i Hl Pl - ( 31 ) 

m j^l 

This last equation admits a minimal solution which is obtained by imposing a 
further condition, namely 

P m Z 1 P n = m= 1,2,... . 
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For n > 1, we will adopt an analogous reasoning. Indeed, given an operator 
X, let us set 

g n (x; z u ..., z n ) ^Y^ 1 —. Y, ad ^ ■ ■ ■ adz ^ *> ( 32 ) 

' — ' to! c — ' 1 

m—1 ki^ \-k m — n 

with n > 1. Then we can define the operator function 

n-l 

Gn(i?o, ■ • • , H n \ Z\, .... Z n -\) := ^ ' Gn-m(H m ; Zi, .... Z n - m ) 

m=o 

- i[Z n ,H a ] + H n n>2. (33) 

At this point, one can show that the sequence of equations generated by for- 
mula (ffifjl is given by 

C l -i[Z l ,H ]=H 1 , [Ci,JIo]=0 
C n -i[Z n ,H ] = G n (Ho,...,H n ;Z l ,...,Z n _ 1 ) , [C n ,flo]=0 n>2 (34) 



In order to write the general solution of this sequence of equations, it will be 
convenient to introduce a shorthand notation; given a linear operator X, we set: 

d*lk E p mXP m := (\X\) {Pm} , (35) 

m 

|>*<k ^X-dXl^^^P.X^^DXd^j, (36) 

Notice that for the superoperators (| • \} Hq and |) • (\jf , which differently from 
the superoperator [|) • (\] H do not depend on the eigenvalues of Ho, we have 
introduced the respective alternative symbols (| • |){ Pm } and |) ■ {\r P \ that will 
be used in the following whenever a certain set of spectral projections {P m } (in 
general, not eigenprojections) of a selfadjoint operator are involved. 
Now, assume that the first n equations have been solved. Then, the operator 
G„+i(-£/o, ■ ■ ■ , H n+ i; Z\, . . . , Z n ) is known explicitly and hence 

C n +i = (|G n+ i(i?o ; ■ ■ ■ j H n+ i; Zi, . . . , Z n )\) Ha , (38) 



[Z n +i,H ] — i |)G n+ i(i?o, • • • , H n+1 ;Zi, . . . , Z n ){\ Ho . (39) 

Again, this last equation determines Z n+ \ up to an arbitrary operator (|^n+i|)ir 
commuting with Hq] in fact, we have: 

Z n +i = {\Z n+ t\) Ho + [|)G n +i(-Ho! * * ■ s H n+1 ;Zi, . . . , Z n )(\] Hg . (40) 



9 



We stress that, in general, the choice of a particular solution for Z n+ i will also 
influence the form of C n +2, %n+2, ■ ■ ■ • 

Thus, we conclude that the sequence of equations defined above admits in- 
finite solutions (even in the case where Hq has a nondegenerate spectrum). 
However, there is a unique minimal solution { M C n , ^ZnjngN which fulfills the 
following additional condition: 

(rZ n \) H =0, n=l,2,.... (41) 



In order to clarify the link of our approach with standard perturbation theory 
for linear operators, let us recall a few facts (see EH)- It is possible to show 
that, under certain technical conditions, there exist positive constants ri,r2, • • ■ 
and a simply connected neighborhood X of zero in C such that, for any A e I 
and m — 1,2, ... , one has that: 

1) the following contour integral on the complex plane 

Pn(A) = ^-jdz (z-tf(A)r 1 (42) 

— where T m is the anticlockwise oriented circle [0, 2ir] 3 9 E m + r m e 10 
around the eigenvalue E m — defines a projection (P m (A) 2 = P m (A)), 
which is an orthogonal projection for A E X n M, with P m (0) = P m , and 
I3Ah P m (A) is an analytic operator- valued function; 

2) the range of the projection P m (A) is an invariant subspace for H(X) (but, 
if the range of P m is not 1-dimensional, in general not an eigenspace), 
hence 

H(X) P m (A) = P m (A) H(X) P m (A); (43) 

3) there exists a (non-unique) analytic family A i— > ^(A) of invertible oper- 
ators such that 

Pm = T^(A)- 1 P m (A) W(X), W(0) = Id (44) 

— with Vt^(A) unitary for A real — which is solution of a Cauchy problem 
of the type 

iW'(X) = J(X)W(X), W(0) = Id, (45) 

where the apex denotes the derivative with respect to the perturbative 
parameter and A i— > J(X) is any analytic family of operators — selfadjoint 
for A real — such that 

D J ( A )(l{p m (A)} =i^P m (A)P m (A) * P m (A)P m (A)P m (A) = 0*. (46) 
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In standard (Rayleigh-Schrodinger-Kato) perturbation theory, one can obtain 
the perturbative corrections to unperturbed eigenvalues and eigenvectors ex- 
ploiting (see, for instance, ref. [201) formula 142(1 and the expansion of the re- 
solvent operator (z — if(A)) , namely 

(z-H - Ho(X))- 1 = {z- Ho)- 1 - X (z - Ho)- 1 Hi (z - Ho)- 1 

- X 2 (jiz-H )- 1 H 2 {z-Ho)- 1 

+ {z- Ho)- 1 ^ (z - Ho)- 1 ^ (z - Ho)' 1 ) +... . (47) 

In our approach we use, instead, properties 2) and 3). Indeed, let us define the 
operator H(X) by 

77(A) := WiX^HiX) W{X), (48) 

which, for real A, is unitarily equivalent to H(X). Using relations 143(1 and ((44(1 . 
we find 

77(A) P m = W(X)- 1 H{X)P m (X)W(X) 

= IU(A)- 1 P m (A) H(X) P m (X) W(X) (49) 

and hence: 

77(A) P m = P m 77(A) P m m = 1, 2, . . . . (50) 

It follows that 

[77(A), H ] =0 (51) 
and then we obtain the following important relation: 

[IU(A)- I i/(A) W(X) - H 0l H Q ] = 0. (52) 

Thus, if we set 

W(X)- 1 H{X)W{X) -H Q = C(A), W{X) =cxp(-iZ(A)), (53) 

and we apply relation 1(13(1 . we find precisely formula ((241) . 

Concluding our treatment of the time-independent case, it is worth stressing 
that, due to conditions 1 (21(1 . for the overall evolution operator we have: 

U(X;t) = e -iHot e -iZ{X;t) e -iC{X)t e z W 

* Z(A; t) = e tHat Z(X) e- lHat * = e -4 *^ e- lHat e - lC{x)t e lZW 

* [C(\),H ] = 0* = e -iZ(*t) e -i(H +C(\))t e iz(\) ^ ^ 

or, more explicitly, 

U(X;t) = e~ iZ ^ ^ e -HE m +cW)tp me izw 

m 

= e~ iZ ^ J2 c M-i(E m + lC\ rn (X))t) P m e iZ ^\ (55) 

m 

where we have introduced the 'reduced rank operators' 

LCJ ro (A) := C(X) P m = P rn C(X) P m , m = 1, 2, . . . . (56) 
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3.2 The adiabatic approximation 



Let us now consider the case where the perturbative adiabatic approximation can 
be applied. This approximation consists essentially in partially neglecting the 
last term in the r.h.s. of equation l|2(J|l — the one involving the time derivative 
of (A, i) i— > Z(X;t) — under suitable conditions. 

Precisely, we will assume that the unperturbed Hamiltonian has (instanta- 
neously) a pure point spectrum with time-independent eigenprojectors. Namely, 
we will suppose that there exists a set of orthogonal projectors {P m } m =i,2,... 
forming a resolution of the identity — Id = J2 m ^ m — sucn that ^ coincides 
with the set of eigenprojectors of H (t), for any t. 

This hypothesis prevents the possibility of occurence of 'level crossings' in the 
spectrum of the unperturbed Hamiltonian. Indeed, it implies that there exist 
real functions 1 1— * Ei(t), t i— ► E 2 (t), . . . such that, for any t, {E m (t)} m= i^,... is 
the set of the eigenvalues of H (t) — specifically: H Q (t) P rn = E m (t) P rn — and 
hence 

Ex(t) + E 2 (t) ? . . . Vt. (57) 



We will further assume that the functions {t i— > £? TO (t)} m =i,2,. 
Moreover, notice that — since [H (t), H (t')] — 0, Vt,i' - 
evolution operator will be given by: 



U (t) = e -*/o t Ho(t)dt = 



i/o E m {t)dt p 

1 jr 



belong toC^R). 
the unperturbed 



(58) 



Then, we will consider a solution {t i— ► C(X;t), 1 i— > Z(X; t)} of the equation 



E ^ a 4(A; t) (^o(i) + ffo(A; t)) = H (t) + C(X; t) 



^ (fc + 1)! "Yo'CCAitJdt 



C(A;t) 



+ 



E 

fc=0 



■ad? 



(fc + l)!~*(Ai*) Z(A;t) 



, (59) 



H (t) 



subject to the condition that 

[C(X;t),P m ] = Vt, Vm 

and 

^ (fcTT)I ad ^)^ (A;0 



fe=0 



(60) 



(61) 



H Q (t) 



We observe explicitly that such a solution may not exist. In more detail, one 
can show that a solution of equation (|59() verifying condition i|60[l alone always 
exists (and, as in the time-independent case, it is not unique), but, in general, 
not a solution verifying condition l|61|l too; we do not insist on this point here 
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since it will be clarified afterwards. We will then say that the perturbative 
adiabatic approximation is applicable if a solution of equation 159fl verifying 
both conditions ((61) (1 and 1(61|) does exist. 

Condition ((6U(I is analogous to condition ((23(1 which has been assumed in the 
time-independent case. Indeed, it is equivalent to the condition that 



[C(X;t),H o (t')} = W, Vt'. 
Integrating with respect to time relation l(6U|) . we find also that 



C(A; t) dt,P„ 



= Vt, Vm. 



(62) 



(63) 



Then, as a consequence of relation (|60|l . recalling formula iJSSJ, we have that 



therefore: 



C(X;t),U (t'] 



U (t),C(X;t) + J2 



C(X;t) dt, U (t') 



Vi,t'; 



- (fc + 1)! io c (V) dl 



C(A;t) 



= Vt. 



(64) 



(65) 



Now, observe that if there exists a solution {f i— > (7(A;t), i i— > Z(A;i)} of 
equation 1(59(1 verifying condition (|60|l (hence relation (165(1 ) and condition (|61() . 
then equation 120(1 will be approximately satisfied setting C(A;t) = (7(A;t) and 
Z(X;t) = Z(X;t). Thus equation 1(59(1 will be the starting point for obtaining 
a perturbative expansion of the evolution operator when the perturbative adia- 
batic approximation is applicable. This task will be pursued in the next section. 
Since at this stage the approximation considered may appear as a mere ad hoc 
computational expedient, we will devote the last part of this section to showing 
what is its meaning; in particular, why we call it 'adiabatic'. 
We will first show that the evolution operator U(X; t) associated with the solu- 
tion {t i— > C(X]t), 1 1— > Z(X;t)} behaves like an adiabatic evolutor. To this aim, 
let us set 

W(X;t) := exp(-iZ(X;tf) , 
and let us define the projection 



P m {X;t) :=W(X;t)P m W(X;t)- 1 
and the selfadjoint operators 

H) fc 



(66) 
(67) 



C(A;t) :=C(X;t)+J2 



ad?**, 



k=l 



(k + l)\ aU /oC(A;t)dt 



C(A;t), 



(68) 



K(X; t) := i W(X; t^WiX; t)=J^ 



adl 



fc=0 



(A + 1)1 z ^ 



Z(X;t) 



= K,(X;t)+jC(X;t), 



(69) 
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with 



£(A;t) 



£(A;i) 



fe=0 v ' 



(fc + 1)! 



ffo(t) 



(70) 



The operator £(A;i) (by virtue of relations l|60|l and l|63|l) and the operator 
/C(A; t) (by definition) satisfy: 

[&(A; i), P m ] = 0, [£(A; f), P m ] = 0, Vi, Vm. (71) 
Then, equation (|59J) can be rewritten as 

H(X; t) = W{\; t) \ H (t) + £(A; i) + K(A; t)J W(A; i)" 1 (72) 

and, by this equation, we have: 

H(X; t) P m (X; t) = W(X; t) \H (t) + £(A; t) + £(A; t)j P m W(X; i)" 1 

* relations |ZU> * = W(X;t)P m \ H {t) + £(A;t) + £(A;t)J P m W(A;t) _1 

= P m (X;t)H(X;t)P m {X;t), (73) 

i.e. i/ie range of the projection P m (X; t) is an invariant subspace for H(X; t). At 
this point, it is convenient to make a detour. 

We remark that as in the time-independent case — instead of assuming that 
decomposition (I72[l holds — one can prove that, under suitable conditions, there 
exists a complete set {P m (A;i))} m =i,2,... of spectral projections of H(X;t) such 
that the functions A i— > P m (A; t) and 1 1— » P m (A; t) are analytic and 

P m (0;0) = P m , 771 = 1,2,... . (74) 

Then, there is a (non-unique) unitary operator W^A; t) such that 

P m (X;t) = W(X;t)P m W{X;t)- 1 , Vi, Vm. (75) 

The operator PF(A;i) can be decomposed as 

W(A;t) = i4(A;t)W(A), (76) 

where A i— > W(A) and (A, i) i— » A(A;i) are unitary operator-valued functions 
such that 



P m (A;0) = W^(A)P m W(A)- 1 , 

P m (A;t) = A(A; t) P m (A; 0) A(A; i), Vi, Vm. 



(77) 
(78) 
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One can show that these functions are solutions of Cauchy problems of the type 
iW'(X) = J(X)W(X), W(Q)=ld, (79) 



i A(A; t) = if (A; t) A(X; t), A(X; 0) = Id, (80) 
with J(A) and K(X;t) selfajoint operators satisfying the following conditions: 

i>^( A )(i { p m( A;o)} = ^ p ™( A ;°) p ™( A ;°) 

m 

* = ix Id = ^A^ Pm(A; * )2 * = -'E^^ )^^ )' ( 81 ) 

m m 

K(X;t) := |)X(A;t)(| {Pm(A;t)} - i^P m (A;i)P m (A;t) 

rn 

* = ^Id= l^p m (A;t) 2 * = -i^P m (A;i)P ro (A;t). (82) 

m m 

Conversely, if A i— > W(X) and (A, £) i— > A(A; £) are operator-valued functions 
which are solutions of the Cauchy problems l|79[l and l|80|l (with J(A) and K(X; t) 
selfadjoint operators subject to conditions l|81l) and (|82[l ). then they will satisfy 
equations (O and J75J|. 

In fact, suppose that relation (|78|1 is satisfied (we will only argue for A(X; t) since 
the argument for W(X) is analogous). Then, setting K(X; t) = i A(X; t) A(X: £) -1 , 
we have: 

iP m (A;t) = [i<r(A;i),P ro (A;t)], Vt, Vm. (83) 
Conversely, if this relation holds, we have: 



A(A; t)- 1 P m (A; t) A(A; t) - A(X; t)- x P m (X- 1) A(X; t) + A(X; t)~ 1 P m (A; t) A(A; i) 

- A(A; i)" 1 A(X; t) A(X; t)' 1 P m (A; i) A(A; t) 
= 0, Vt, Vm, (84) 

where we have used the fact that A{X]t)^ 1 = — A(X; t)~ 1 A(X; t); hence: 

A(X; i)- x P m (A; t) A(X; t) = A(A; 0)- 1 P m (A; 0) A(X; 0) = P m (A; 0), Vt, Vm. (85) 

Thus, relation (|78|l holds if and only if K(X;t) satisfies equation JS3J). At this 
point, using the fact that P m (A; t) P m (A; t) P m (X; t) = 0, Vm, one can check 
easily that 

K(X;t) = K(X;t) + K(X;t), (86) 

- with K{X-t) := {\K{X-t)\) H(x . t) and #(A;t) := \)K{X;t)(\ H(x . t) - solves 
equation (|83|l if and only if: 

K(X;t)=iJ2Pm(X]t)P m (X;t). (87) 
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In conclusion, it is proven that, if there exists a complete set of spectral pro- 
jections {P m (A;t)} m= i ) 2,... of H(X;t) with the properties specified above, then 
equation l|75|) is satisfied by only and all the unitary operator-valued functions 
(X,t) i-> W(X;t) such that W(X) = W(X;0) and A(X; t) = W(X)~ 1 W(X; t) are 
solutions of the Cauchy problems (|79fl and (JSDJ, with conditions (I81|l and (|82|) . 
Moreover, if one sets 

H{X: t) := W(X; t) _1 i/(A; t) W(X] t), (88) 
the following relation holds: 

H( A; t) P m = P m H(X; t) P m , Vi, Vm. (89) 



Then, setting 



K(A;t) = i( iy(A;t)" 1 T^(A;i) ) , (90) 

\ / H(j(t) 

£(X;t) =H(X;t)-H (t)-!C(X;t), (91) 



one finds that 



[£(X;t),P m ]=0, Vi, Vm. (92) 

In the next section, it will be shown that equation (|68fl allows to recover the 
operator C(X;t) from C(A;i). Thus one can actually define C(X;t) through 
formula (|68|l . One can also prove that if relation l|92|l holds then, as a conse- 
quence, [C(X; i), P m ] = 0, Vt, Vm. Hence, setting W(A;i) = exp(— i Z(X; t)), 
one reobtains equation (|59(l with condition H6Q(I automatically satisfied. 

Let us now come back to our original purpose of investigating the behaviour 
of the evolution operator U(X;t) that approximates U(X;t): 

U{X;t) = Uo ( t ) e -izW e -imx;t)dt e iz W 

* Z(X; t) = Ad Uo(t y Z{X; t) * = e ~ lZ{x ^ U (t) e~ l ^ c ^) dt e lZW 

* adiabatic approximation * ~ e - lZ(X:t) U (t) e'^ 6 ^ dt e tZ(x) 

= :U(X;t). (93) 

Recalling formula l|58|) and relation the adiabatic evolutor U(X;t) can be 
written as 

U(X;t) = e^ z(x ^ ^ e -/o t (s m (t)+c(A;t))dtp me iz(A) ! (Q4) 

m 

or — introducing the reduced rank operator 

[C\ m (X;t):=C(X;t)P m =P m C(X;t)P m , m=l,2,... (95) 
— in the more expressive form: 

U(X;t) = e -^(^t) ^ e -^ t (B m (t)+LCJ m (A;t))dt Pmei 2(A)_ (Q6) 



1G 



To show that U (A; t) behaves indeed as an adiabatic evolutor, let us observe 
that 

C/(A;i)Pj(A;0) = e ~ a ^ U (t) e^o^) dt e iJ W Pj(X; 0) 

nt 

= e -iZ(\;t) p e -i/ t (£; ; ,(t)+C(A;t))dt el i(A) 
= e -^(A;t) p. Uq ^ e - l / t C(A;t)dt e ^(A) 

* definition (S3* = Pj (A; t)U(X;t), j = l,2,...; (97) 

namely, U(X;t), as it should, intertwines the projection P,(A;0) with the pro- 
jection Pj(X;t). Let us do the following observations: 

• The evolutor U(X;t) intertwines spectral projections, that in general, for 
A / 0, are not eigenprojections. This situation is more general than the 
one considered originally by Kato [23 an d earlier by Born and Fock in 
their seminal paper |22j . Nevertheless, due to its importance in several 
applications, this situation has been studied in later times by other authors 
(see, for instance, Nenciu's paper [231 )■ Considering this more general 
situation is in our case unavoidable, since in presence of the perturbation 
(i.e. for A ^ 0) the unperturbed eigenvalues ('energy levels') can 'split'. 

• For A = 0, the adiabatic evolutor reduces to the unperturbed evolution 
operator: 

U(0;t) = U (t). (98) 
This fact justifies the term ' perturbative adiabatic approximation'. 

Anyway, in order to better highlight the typical structure of an adiabatic 
evolutor, it is convenient to rewrite the expression of U (A; t) as follows: 

U(X;t) = A(X;t) J2e~ iS ° {Emit)+Um(X ' ,t))dt Prn(*0), (99) 

m 

where 

A(X; t) := W(X; t) W(X; 0)" 1 , 

n m (X; t) := W(X; 0) C(X; t) W(X; 0)" 1 P m (A; 0) (100) 

= W(X;0) LCJ m (A;t)^(A;0)- 1 . 

Assume that the time-dependence of the Hamiltonian H(X; t) is characterized 
by a time scale T > 0, i.e. 

H (t) = Mo(t/T), H«(X;t) =Mo(X;t/T), <e[0,T], (101) 

with [0, 1] 3 s i ► JGq{s) and [0, 1] 3 s i— > ^E (A; s) given operator- valued func- 
tions of the scaled time s. From the physicist's point of view, the parameter T 
measures the 'slowness' with which the non-isolated quantum system described 
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by the Hamiltonian H(X;t) is influenced by the external world. Then, for the 
spectral projections {P m (A; i)} m =i,2,... of the Hamiltonian H(X;t) we have: 

P m {X;t) = V m {X;t/T), Vi S [0,T], Vm, (102) 

where 'P m (A; s) is a projection- valued function of the scaled time. 

Recall now that the operator A(X; t) is solution of an initial value problem of 

the form 

i A(X; t) = K(X; t) A(X; t), A(A;0)=Id, (103) 
where K(X;t) has to satisfy the following condition: 

IW;*)<l{P m (A !t » : = *E Pm(X;t) P m (X;t) 

m 

* V m (X; s) = d s V m (X; s) * = i^P m (A;i/T)7> m (A;t/T). (104) 

m 

Besides, as K(X;t) = i W{X; t) W(X; recalling definition we find that 
the selfadjoint operators K{X\t) and IC(X;t) are unitarily equivalent since they 
are linked by the following formula: 

JC(A; t) = W(X; t) _1 i^(A; t) W(X; t). (105) 

More specifically, using the fact that 

\) Ad w(\:t)-i(-)(\ Ho(t) = |) Ad TV(A;t)-i(-)(| {Pm} = Ad W ( A;t) -i |) ■ (| {Pm ( A ;t)} , (106) 
we have: 

£(A; t) = W{X; t) _1 ^(A; t) W(A; t) = O f ij . (107) 

This observation 'completes the picture'. Indeed, it turns out that, in the adi- 

abatic limit, the contribution of the operator K.(X;t) in equation (|20|l can be 
neglected, namely that approximation (|61|l is justified as it is equivalent to the 
standard adiabatic approximation. 



4 Solution with the perturbative adiabatic ap- 
proximation 

Let us now face the task of providing the perturbative solutions of equation (|59|) , 
which we recall is subject to condition (|60|l and to what we can at this point 
legitimately call 'adiabatic approximation', i.e. condition (|61ll . 
It will be convenient to express equation Q59J1 in terms of the operator £(A; t) 
in place of C(X;t). In fact, as it will be shown later on in this section, one 
can develop a simple perturbative procedure which allows to compute, order 
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by order, the operator C(X;t) from <t(\;t). Then, the equation that has to be 
solved perturbatively is the following: 



oo .u 



E J\ a 4(A; t) W) + ^o(A; t)) = £(A; t) - ff<>(A; i) 
fe=i 



E 



ad"; 



+ 1)! Z ( A - 4 ) 



,2(A;i) 



,(108) 



where the operator <£(A; i) is subject to the costraint 

[€(X;t),H o (t)}=0, W. (109) 

It will be also shown that condition l|109(l is actually equivalent to condition (ffjQfl 
of which it is a straightforward consequence. This fact is not immediately ev- 
ident from the definition of <£(A;i). We will postpone the problem of checking 
the validity of the adiabatic approximation (condition JBTJ) to later analysis. 
Given linear operators X, X\ , . . . , X n and Y\ , . . . , Y n , let us set 

(±i) m 



T^n • • • i 



E 



(m + 1)! z - 

m=l v ' fciH hfc m =" 

Then, for n > 2, we can define the operator function 



]T ad y&i • • • ad nm X, n>l. (110) 



n-l 



R« ■ ■ ■ 7 x n -i, Yi, . . . , i^-i) :— ± iz^_ rn (x m ;Yi, . . . , y„- to ). (in) 

Noting the analogy of equation ll 108(1 with equation (|24|l obtained in the time- 
independent case, we conclude that, given the power expansions 

oc oo oc 

H4X;t) = J2 xnH n{t), £(\;t) = ^\ n £ n (t), Z(\;t) = J2 xn Zn(t), (112) 

n— 1 n— 1 n—1 

the sequence of equations generated by formula l|108|) . with condition (|109|l . has 
the following form: 



Ci (t) - i ii (t) , ff (t) = J?i (t) - ( Zi (t) 



ffo(t) 



0, 



(113) 



£»(<)-» Z n (t),H (t) 



= G„(H (t), . . . , ff n (t); . . . , £„_i(t)) 

R+ (i x (t ) , . . . , L- 1 (t ) ; Z x (t ) , . . . , Z n - x (t ) ) 

z n {t) 



£„(*), H (t) 



= 



H (t) 

n > 2, 



Ho(t) 



(114) 
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Exactly as in the time-independent case, this system of equations can be solved 
recursively. Indeed, we have that 



C X (t) = (\Hl(t)\)H (t)- u(\Mt)\)H (t), 

Z x {t) = (\Zi(t)\)H (t)+l\)Hi(t)(\} Ho (t), 



(115) 
(116) 



— where (\Zi(t)\)H (t) is an arbitrary operator commuting with Ho(t) (notice 

that ^(l-^i (*)l)-ffo(t) = (l-^i Wl)-ffo(t)) — &nd, assuming that the first n equations 
have been solved, so that the operator functions 

G„+i(. . . , H n+1 (i); ...,Z n (t)) and R+ +1 (. . . ,i n (i); . . . , Z n {t)) 
are known explicitly, the solution of the (n+l)-th equation is given by 

£n+i(t) = ( | G„+i (Ho (t) , . . . , H n+ i (f ) ; Z x (t) , . . . , Z n {€)) 
R++i(ii(i), • • • , L(t); Zi(t), . . . , Z n (t)) ) 

4+i(t)L,. ( 117 ) 



d 

dt 



ffo(*) 



•z„+i(t),flo(*) =* )G n+ iCff (*),■■■, ffn+iWsii (*),..., i„(t)) 



(118) 



Again, equations Ijll7|l and (|118|l determine £„+i(i) and Z n +i{t) up to an op- 
erator (|Z„+i(t)|)^f ( t ) commuting with Ho{t) and we have: 



Z n+1 (t) = ( Z, i+1 (i) 



H (t) 



G n+1 (H (t),..., H n+1 (t) ; Zi (t) , . . . , Z„ (i) ) 



H (t) 



(119) 



Eventually, one has to check that condition 1)61(1 is satisfied. Explicitly, one has 
to check that the solution obtained {t i— > C„(A; t), 1 1— > Z„(A; t)}„eN is such that 



0. 



R+(Zi(t), . . . , Zn-iWiZttt), . . . , Z n ^(t)) - Z n (t) 



ffa(i) 



0, n > 2, (120) 



At this point, we have to show how the operator C(X;t) can be recovered, 
at each perturbative order, from the operator £(A; t). To this aim, let us recall 
that 

C(X; t) = C(A; t) - £ ad^ (A;t) dt C(A; i). (121) 
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Now, notice that, if we substitute in equation i|121l) the power expansions 

C(X; t) = T,n=i Xn Cn(t) and t(X;t) = J2n=i A ™ £«W > and we sin S le out thc 
various perturbative orders, we conclude that the n-th order, which on the l.h.s. 
is given simply by A™ C n (t), consists on the r.h.s. of X n £„(i) plus a function of 
Cx(t), . . . , C n (t) and J Q Ci(t) dt, . . . ,J Q C n -i{t) dt. Thus, exploiting this fact, 
we can achieve an order by order solution. Indeed, recalling definition 
one finds out that C(X; t) can be obtained order by order from <£(A; t) by means 
of the following recursive process: 

Ci(t) = Ci(t), 

C n (t) = R-(d 1 {t),...,C n - 1 (t);^ Cx(t) dt,...,^ #„-i(t) dt 

+ l n (t), n > 2, (122) 



This recursive process allows to easily prove by induction that the commutation 
relation [£(A; t), H {t)] = (or, equivalently, [£(A; t) , P m ] — 0, Vm) implies 
that [C(X;t),H (t)} = 0. Thus equation (|108|) . with condition (|109|l . is indeed 
equivalent to equation (|59|l . with condition l|60|l . 

We want to show now that there is also another recursive process allowing to 
recover C(A; t) from £(A; t) which less expensive from the computational point 
of view. 

To this aim, let us define the function avxp : M — > R + in the following way: 

I r x e x - i 

avxp(x) = - / e* dt = for x ^ 0, avxp(0) = 1. (123) 

x Jq x 

This function extends to an entire holomorphic function on the complex plane 
which is given by 

avx ?W = E TfcTTv A (124) 

fc=0 v *' 

We recall that the function avxp is fundamental in the theory of Lie groups 
since it is strictily related to the differential of the exponential map, usually 
denoted by dexp (see ref. |H]). ^ n fact, given a Lie group G with Lie algebra g, 
identifying the tangent spaces at any point of G and of g with g itself, one has: 

dexppf)y = avxp(-X)Y", X,Y e g. (125) 

Besides, it is well known (see, for instance, ref. [5S]) that the the meromorphic 
function 1 /avxp admits the following expansion in the open disk of radius 2tt 
centered at zero: 

avxp(z)- 1 =£f **, (126) 

fe=0 
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where {/3q, fli, ■ ■ .} are the Bernoulli numbers, namely the rational numbers 
defined recursively by 

A) = l, ( k t 1 )l3o + ---+ ( h t ^ Pk = fc = l,2,.... (127) 



o y ^ u 1 v k 

We recall also that 

/W=0, ^ = (-l) fc+1 , fc=l,2,.... (128) 

|P2fe| 

Now, notice that, according to definition l|68(l and formula (|124[l . we have: 

C(A; t) = avxp(-i adj*^ dt ) (A; f). (129) 
Then, by means of formula l|126fl . we can write 

C(X;t) = avxp(-iadj ot( 5 (A . l)dt ) £(A;i) 

= g(z^ ad ^ (A;t)dt £ (A . t) . (130 ) 

Let us observe that, as in the previous case, if we substitute in this equation the 
power expansions C(X;t) — Yl^Li A™ Cn(*) an d <£(A;t) = 2~Z^Li A" &n(t), and 
we single out the various perturbative orders, we find that the n-th order, which 
on the l.h.s. is given simply by \ n C n (t), consists on the r.h.s. of a function of 

<£i(t), . . . , <t n {t) and J Q Ci{t) dt,...,J Q C n -i(t) dt. Thus, again, we can achieve 
an order by order solution. 

To this aim, given linear operators X and Yi, . . . ,Y n , n > 1, let us set 

Bt(X;Y 1 ,...,Y n ):= T Y, L E ad rfci • • • ad^ m X; (131) 



namely, we have that Bf(X;Yi) — ^ady 1 X (since Pi = —1/2) and, by virtue 
of relations l|128|) , 

Bt(X;Y u . . . , Y n ) = l - & d Yn X± lJ ^T E ad ^ ' ' ' ad ^ X ' 

m=l ' k!-i \-k2m=n 

(132) 

for n > 2, where p(n) is equal to n if n is even and to n — 1 otherwise. Next, 
for ?i > 2, given linear operators Xi, . . . , A n _i and Yi, . . . , Y n -i, we can define 

ri-1 

B^(Xi, . . . , X„_i; Yi, . . . , Yn-i) := ^ B^_ m (X m ; Yi, . . . , Y„_ m ), (133) 

m— 1 
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where the definition of B+(Xl, . . . , X n _x;Yi, . . . , will be used in the next 

section. Then, one finds out that C(X; t) can be obtained order by order from 
£(A; t) by means of the following recursive process: 

<?i(i) = &i(t), 

d n (t) = B-(t 1 (t),...,c r ^i(t);J (?i(t) dt, . . . , ^ (? n _i(t) dt 

+ C„(t), n > 2, (134) 



One can verify that this recursive process is computationally convenient with 
respect to the recursive process l|122[) . Specifically it turns out that, from the 
4-th perturbative order on, each step involves a smaller number of terms in 
comparison with process (I122[l . with a gain which increases step after step. 



5 The general case 

In this section we will consider equation l|15l) in its full generality. This equation, 
taking into account definition JBSJ, can be rewritten as 

00 -k / i \ 

E h ad ^;') *) - feTT i(A; t] = £(A; tl (135) 

k=0 ' ^ ' 

where the operator C(A; t) is linked to the operator C(A;i), which appears in 
decomposition JSJ of the interaction picture evolution operator, by the same 
relation that links the operators £(A; t) and C(X; i) introduced for studying the 
adiabatic approximation; we rewrite it here for the sake of clarity: 

£(A;t) := avxp^-i^ C(A;t)dt) C(A;i) 

= C ^ *) + E TT^TV 4c ( a ;1) dt C(A; t). (136) 
fc=i v 

We already know (see section[3J) that the operator C(A; t) can be recovered from 
the operator £(A; t) by means of an order by order procedures. Thus, we can 
solve equation l|135f) for £(A; t) up to a given perturbative order and obtain the 
perturbative expansion of C(X;t) truncated at the same order parallely. 
First of all, we want to provide an interpretation of decomposition JHJ which 
sheds light on its meaning. To this scope, notice that this decomposition can 
be rewritten as 

T C (A; t) = T Z {X- i)t T(A; t) T Z (X; 0), (137) 
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where 

Tc(X;t) := exp^-i^ C(A;t)dtj and T Z (A; t) := exp(-i Z(A; t)) . (138) 

Now, formula (|137|l can be regarded as a passage to a further (generalized) 
interaction picture performed on the Hamiltonian H(X;t). Indeed, assuming 
that Tz(X;t) satisfies the equation 

if z {X;t) =T z (A;t)3(A;t), (139) 

one finds that 

* T C (A; t) - (t z (X; i)t tf(A; i) T Z (A; t) - 3(A; t)) Tc(A; i)- (140) 
Then, since 

3( A ^) = E7fcTTT! ad W)^(A;i), (ui) 

equation l|135|) expresses precisely the fact that £(A; t) is the transformed Hamil- 
tonian obtained by switching to this new interaction picure; namely: 

C(A; i) = T Z (A; f)t H(X; t) T Z (X; t) - 3(A; t). (142) 

It follows that 

" ' 1 (143) 



T c (X-,t) =exp(-if^X n ^ C n (t)dtj 



is nothing but the Magnus expansion of evolution operator associated with the 
new interaction picture Hamiltonian £(A;t). 

Let us now investigate the perturbative solutions of equation (|135J) . Substi- 
tuting the power expansions 

OO OO OO 

H(X;t) = J2^ n H n (t), <£(\',t) = jr i \ n <t n (t), Z(X;t) = J2* nz n{t), (144) 

n— 1 n—1 n—1 

in equation (|135|L we obtain an infinite set of coupled equations that allows to 
compute the operators {£ n (i)}n6N, {Z n (t)} n &i- In fact, in analogy with the 
time- independent case (compare with definition given linear operators X 

and Y, for n > 1, let us set 

G n (X,Y;Z 1 ,...,Z n ):=J2—, E ad z • • ■ ad z X — . (145) 

* — ' to! ^-^ 1 m \ TO + 1 / 

m=l feiH hfcm=n 



Then we can define G„(#i(t), . . . , #„(*); Zi(i), . . . , Z„_i(t); Zi(t), . . . , £„_i(t)) 

n-l 

53 g n - m {H m {t),Z m {t);Z x {t\... 

i — rn 

{t)) + H n (t), n>2. (146) 



as 

n—1 



m— 1 
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With these notations, one can write the sequence of coupled equations which 
gives a perturbative solution of equation i|135|) as follows: 

Zi(t) =#i(t)-e:i(A;i), 



Z n (t) = G„(#i(i), . . . , # n (t); Zi(i), . . . , Z n _i(t);Zi(t), . . . , ^„_i(t)) 

- n > 2, (147) 



It is clear that, as in the time-independent case, this infinite set of equations 
can be solved recursively. Moreover, recalling the recursive process j 13411 that 
allows to obtain, at each perturbative order, the expression of C(X;t), one can 
calculate order by order both the operators {C„(t)} ne N and {Z n (t)} n( z^. In- 
deed, integrating with respect to time each equation in the sequence (114711 and 
combining the new sequence of equations so obtained with the recursive pro- 
cess ((T34) (or (H220, we find: 

rt 



Zi(t) = J (£i(t)-ei(t))dt+Zi, 

Ci(t) = Ci(t), 



Z„(t) = J (e w (... > fl B (t);...,Z B _i(t);... > Z B _i(t)) -<£„(*)) dt+Z„, 
C„(t) = B-^...,£ n - 1 {t);...,J C»_i(t) dt) + £»(*), n>2, (148) 



This time, differently from the equations obtained in section [3] for the time- 
independent case, at each perturbative order we have a couple of equations. 
The solution of the first couple of equations is obtained by simply choosing the 
arbitrary operator- valued function t i— > £i(i) and the arbitrary operator Z±; 
similarly, the solution of the n-th couple of equations, for n > 2, involves the 
previously computed functions 1 1— > Z\(t), . . . , i i— » Z n —\(t) and requires only the 
choice of the arbitrary operator- valued function i i— ► £„(i) and of the arbitrary 
operator Z„. This choice can be fitted according to computational convenience. 
Notice that, in particular, the choice of the operators {Z n } n ^ determines the 
initial condition 

oo 

Z(A) = Z(A;0) = ^A"Z n . (149) 

We will see soon that one can set quite natural conditions which fix a unique 
choice of these arbitrary quantities uniquely. Before doing this, we want to 
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show that it is possible to write a sequence of equations which is equivalent and 
structurally similar to the one given above but such that the solution of the 
n-th couple of equations — the one determining C n (t) and Z n (t) — does not 
involve explicitly the operators Z%(t), .. . , Z n -i(t). 
The first step is to rewrite equation l|135|l as 

avxp(i ad z(A;4) ) Z(X; t) = exp(i ad z(A;4) ) H(X; t) - C{\). (150) 

Next, using the relation 

avxp(z) -1 exp(z) = avxp(— z) _1 , (151) 

we have: 

Z(X; t) = avxp(i ad Z (x-t)) 1 (cxp(i ad z(A;t) ) H(X; t) - C(X)j 

k=0 

Recalling the fact that (3q = 1, Pi = — 1/2 and relations (|128fl . we can further 
simplify equation (|152|l : 

Z(X; t) = H(X; t) - C(X) + ~ ad z(A;t) (#(A; t) + C(A)) 

+ E(-!) fc ^ ad ^(A;*) (H{X- 1) - C(A)) = | ad z(A;t) (#(A; t) + C(A)) 
fc=i 

+ ( x - E ^ ad2 4^)) (^( A ; *) - C ( A )) • ( 153 ) 

Notice that in the time-independent case (H (X; t) = H (X)), assuming as in 
section that [C(X),Hq] = and Z(X;t) = exp(iHot) Z(X) exp(-iHot), hence 

Z(X; t) = -i Ad cxp{lHot) ad z(A) H Q , (154) 

from equation (|153fl . using relation (|18|l . we find: 

- i ad z(A) H = ffo(A) - C(X) + l - ad z(A;t) (ff (A) + C{X)) 

- E w ad ^*) - • ^ 155 ) 

fc=i 

This equation yields another order by order solution procedure with respect to 
the one described in section Anyway, we will not insist on this point and we 
leave the details to the reader. 

Returning to the general time-dependent case, from equation l|153|) . we can 
obtain an infinite set of equations which allows to give perturbative expressions 
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of Z(A; t), C(A) and can be solved recursively. To this aim, given linear operators 
X , Y and Z\ , . . . , Z n , let us set 

B n (X, Y;Z U ..., Z n ) := £ £ ad Zfci • • • a d Zk J(-l) m X - Y), (156) 

m=l kt-i \-k m =n 

forn> 1; namely, B\(X, Y\ Z\) = \aAz x {X -\-Y) and, recalling definition 1(132(1 . 

B n (X,Y;Z 1 ,...,Z n ) := l -ad z jX + Y) 
p(")/2 . R . 

+ E E ad, fci ...ad^(y-X) 

"1=1 fciH hbm^n 

= Zi, . . . , Z n ) + B+{Y- Zi, . . . , Z„), (157) 

for n > 2, where p(n) is equal to ri if n is even and to n — 1 otherwise. Then we 
can define, for n > 2, 

B„ = B„(#i(i), H n -i(t); d(t), Zi(*), . . . , Z n _i(t)) 

as 

B„ := J] B 

n — rn 

(H m (t), £ m (t); Zi(t), . . • , Zn-mf*)) 

m— 1 

= B-(. . . , H„_i(i); • • • , ^n-i(t)) + B+(. . . , C„-i(t); . . . , Z„-i(t)). (158) 

At this point, inserting the power expansions ()144|l in equation 1)152(1 and equat- 
ing the terms of the same order in the perturbative parameter A, we obtain the 
following sequence of equations: 

Z t (t) =#!(*)- Ci (A; t), 
d(t) = Ci(t), 

Z„(t) = B B (£Ti(t), .. . , #„_i(t); Ci(t), . . . , C„-i(t); . .. , Z n -i(t)) 
+ H n (t)-€ n (t), 

C n {t) = B r ; fci(t), . . . , C„-i(t); / Ci(t) dt, . . . , / C n _i(i) dt^ 

+ £„(*), n > 2, (159) 



Again, this sequence of equations can be solved recursively and the solution of 
the n-th couple of equations requires simply the choice of an arbitrary operator- 
valued function (t h- > £„(£)) and of an arbitrary operator (Z n ). 
An important class of solutions is determined by the condition 

<£ 1 (t) = £ 1 (0) = £ u ...,£ n (t)=£ n (0) = £ n ,... Vi. (160) 
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This condition is equivalent to the following: 

C 1 (t)=C 1 (0) = C 1) ... ) C n (t) = C n (Q) = C n) ... Vt. (161) 
Moreover, if this condition holds, we have: 

d =£1,..., C n = £„,... . (162) 
Then the solution of the first equation, namely 

Z 1 ({C 1 ,Z 1 };t)= [ ^(tJdt-tCi + ^i, (163) 
Jo 

is fixed by the choice of the 'arbitrary constants' C\ and Z\. Inductively, the 
solution of the n-th equation is obtained by substituting the previously chosen 
arbitrary constants C\, . . . , C n _i and the solutions 

t^Z 1 ({C l ,Z 1 },t),...,t^Z n _ 1 ({C k ,Z k }^;t) 

of the first n — 1 equations — that are fixed by the choice of the additional 
arbitrary constants Z\, . . . , Z n -\ — in the formula 

Z n {{C kl Z k }l =l -t) = J B„(..., J ff„_ 1 (t);...,a i -i;...,Z„_ 1 ({C fc ,Z fc }^Z 1 1 ;t))dt 

H n (t)dt-tC n + Z n , (164) 



which involves the n-th order arbitrary constants C n and Z n . 

Now, as anticipated, we will give suitable conditions that fix the arbitrary 
constants {C„, Z n } n t£N — hence, a solution of the system (I159JI — uniquely up 
to a certain perturbative order N S N. 

To this aim, it will be convenient to introduce the following notations. First, 
for the sake of brevity, let us define 

J„(...) :=B n (...,H n _ 1 {t);...,C n - 1 ;...,Z n _ 1 (t)) + H n (t), n>2. (165) 

Next, given an analytic function A ^ f(X) = X)^Lo ^ n /ni we wm se t 



/ [N] (A):=^A"/„. (166) 

n=0 

Moreover, given another analytic function A 1— > h(X), we will set: 

/(A)«MA) 44 / [N1 (A)=/ l[N1 (A). (167) 

Finally, given t G (0, 00], we will say that a function t 1— » has ^ero average 
over the time span [0, t] if 



1 /■* 

lim - / F(t) dt = 0, 
t-,t t J 



(F(-)}f := hm - / F(t)dt=0, (168) 
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where, obviously, the limit is essential in the previous definition only in the case 

where i = oo. For t < oo, equation l|168l) expresses the fact that the function 

1 1— > F(t) has a 'purely oscillatory behavior' in the interval [0,i]. 

Then, fixed a certain perturbative order N and r £ ]0 , oo[ , we set the following 

conditions: 

CI the following relation holds: 

Z n {r) = Z n , n = l,...,N; (169) 

C2 the operator- valued function t i— > Z n (X; t), n = 1, . . . , N, has zero average 
over the time span [0, r]. 

Condition CI is equivalent to the following: 

Z im (X;r) = Z im (X). (170) 

Thus, since 

T(A;r) « exp(-iZ w (A;r)) exp(-z C m (A) r) exp(zZ [N] (A)) , (171) 
condition CI implies that 

x N 

T(X;t) rs exp(-iZ [H] (X)) exp(-z C [N] (A) r) exp(i Z [N] (A)) , 

w exp(-i Z(A)) cxp (-i C(A) r) exp(i Z(A)) . (172) 

If we are able to show that condition CI can indeed be satisfied, this result has 
the following interpretation. There exists a hermitian operator $j(X), depending 
analytically on the perturbative parameter A, such that the 1-parameter group 
of unitary operators generated by it interpolates, up to the N-th perturbative 
order, the evolutor T(A; t) at t = r; namely: 

T(A;r) w exp(-i£(A) t) . (173) 

Indeed, this relation is satisfied if we set 

55(A) = expH Z{X)) C(A) exp(i Z{X)) . (174) 

Now, observe that, denoted by { r C n , T Z n }^ l=1 an item of the first N arbitrary 
operator constants satisfying conditions CI and C2, applying condition CI to 
formula fTMjl yields T Ci = (H(-)) T and 

T C n = — J J„ ^. . . , H n (t); . . . , T C„_i; . . . , Z n -i({ T Ck, T ZkYl = \'i t)j di, (175) 

for n = 2, ...,N. Thus condition CI determines the operators { T C„}^ =1 
uniquely for a fixed N-tuple { T Z n }^ =1 . Suppose that 1 1— > H(X;t) is periodic up 
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to the N-th perturbative order, with period T — H [fl] (X;t) = H [N] (X;t + T), Vt 
— or equivalently 

H n {t) = H n (t + T), Vt, Vn G {1, . . . , N}. (176) 

Then the functions t <— > Zi(Ci, Z x ;t), . . . ,t <— * Z^({Ck, Zk} k Zl;t) are periodic, 
with period T, if and only if C x = T C\, . . . , Cm = t Cn, with r = moT, for some 
mo G N. In fact, we have: 

Z x (C x ,Z x ;t + T) = / Hi(t)dt+ / fli(t)dt- (< + T)Ci + Zi 
Jo Jt 

*H X (-) periodic * = J H x (t) di - t C x + Z x + (J H 1 (t) dt - T C x V (177) 

Hence, the function 1 i— > Z x (C x , Z x ;t) is periodic if and only if 

d = (Hi(O)t = T Ci, for r = m T, m G N. (178) 
Next, reasoning by induction and using the fact that the function 

1 1— > J„ f H x (t), . . . , H n (t); . . . , C„_i; Z X (C X , Z x ;t),..., Z„_i({Cfc, Zk}^ = l ; i)^ 
is periodic, with period T, if the functions 

iii(-), . . . , H n {-)\ Z X {C X , Z x ; •),... , Z„_i({Cfc, Zk}frZl; •) 
are periodic, with period T, one proves the claim. 

It follows that, if t i— > H(X;t) is periodic up to the N-th order, then, setting 
r = moT in condition CI for some nonzero positive integer mo, we have: 

T(A;mT) « cxp(-if)(A)mT), to =1,2,.... (179) 

Condition C2 is equivalent to the condition that the operator-valued func- 
tion t i ► Z m (\;t) has zero average over the time span [0, r]. Besides, applying 
condition C2, namely 

-/ ZJt)dt = 0, n = l,...,H, (180) 
to formula fifity yields T Zi = -(/ ( ) H x (i) dt) T + |r T Ci and 
X = -J|^J„(..., #„(t); . . . , r C„_ i; . . . , Z n _i({ T C fc , T Z fc }j£j; t)) dt) d< 
+ ^r T C„, n = 2,...,N. (181) 

Hence, as claimed before, there is a unique N-tuple { T C n , T Z n }^ =x , determined 
by formulae (1 175(1 and (|181fl . satisfying conditions CI and C2. 
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Now, let us suppose that the following limits exist: 



°Ci := lim T C\ = lim - / HAt) dt, 

T — >OC T — >OC T Jq 



°Zi : = lim 



wo 



#i(t) di)dt+-T°°Ci 



°C N := lim ( - / J N 



j^Cn-i! ■ • • , •Z'N-i({ °Cfe, °°Zk}k = i': t)) dt 
™ (~ \ JXJo^ { - Z N-i({°°C fe , "Zfc}^ 1 ; t)) dt) dt 

(182) 



Of 



Then, the following relations hold: 



1 



lim -Z n ({C k ,Z k }l =1 ;t) = 0, n = l,...,N, 

t— tea t 



lim 1 f Z n {{C kl Z k } n k=1 ;l) dt = 0, 



n = 1, 



N, 



(183) 



(184) 



where Ci = °°Ci, Zi = . . . , C N = °°C N , Z N = °°Z N . Conversely, if rela- 
tions (|183|l and 1|184|) hold for some operators {C kl Z k }^ =1 , then the limits (|182[) 
exist and C\ = °°C\, Z\ = °°Z\, . . . , Cm = °°C|\|, Z^ = °°Z^. Thus we can in- 
clude the set {°°C k , °°-Zfc}fc = i among the sets of operator constants determined 
by conditions CI and C2 if we allow t = oo and we rewrite condition CI: 



CI the following relation holds: 



lim — Z n (t) = lim - Z n , 

t~*T t t—>T t 



(185) 



while condition C2 remains unchanged. Indeed, relation l|185|l reduces to equa- 
tion (jl69Jl if r < oo and to relation (|183|) if r = oo. Moreover, relation ()184ll 
expresses the fact that the function 1 1— > Z n (X; t), n = 1, . . . , N, has zero average 
over the time span [0, oo]. 

At this point, considering the time-independent case — Ho(t) = Hq and 
H <> (X;t) = H (X) — it is natural to ask what is the relation between the so- 
lution associated with the arbitrary constants {°°C n , °°Z n } n ^ (if they exist) 
discussed in this section and the solutions obtained in section As we have 
done in that section, we will assume that the unperturbed Hamiltonian Hq has a 
pure point spectrum E\, E2, ■ ■ ■ and we will denote by P\, P2, ■ ■ ■ the associated 
eigenprojectors. Then, we want to prove that: 

1. for any N 6 N, the limits (|182fl — thus the set of operator constants 
{°°C n , °°Z n } n& N — exist; 
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2. the minimal solution { N C„, M Z„}„ e N of the sequence of equations i|34|). 
i.e. the solution obtained imposing condition JHJ, satisfies the relation 

*C„ = °°C n , M Z„ = °°Z„, VneN; (186) 

3. the operator-valued function t Z(X;t) = YlnLi ^" ^n(*) which verfies 
conditions CI and C2, with r = oo, is such that 

Z n (t) = Z n ({™C k , °°Z k }^l;t) = e iHot °°Z n Vn G N. (187) 

In fact, wc know that {^C n (t) = M C*„, ™Z n (t) = e - lHot ™Z n e lHnt } neN is a 
solution of the sequence of equations (|159|1 with condition (|16H . Observe that 

lim -">%,(<) =0, VneN, (188) 

thus condition CI, with r = oo and N arbitrary, is satisfied. Moreover, given a 
linear operator X, we have: 

X(t) := e iHot Xe- iHot 
*since[(|X|) ffo)J ff ]=0* = <|X|) ffo + e iff °* l)*^ e"^* 

= <m) ff0 + E ei(Ei " Ei)t ^' x ^; ( 189 ) 

hence: 

(X(-)) 00 = (\X\) Ho . (190) 

By condition |(1T)|. it follows that ( M Z n (-))oo = 0, n = 1,2,..., thus condi- 
tion C2, with t = oo and N arbitrary, is satisfied. Then, by the uniqueness of 
the solution satisfying conditions CI and C2, for a given r G]0, oo] and up to a 
certain perturbative order N, we must conclude that 

™C n = °°C n and ^Z„(<) = Z„({ 00 C fe , 00 Z fc }^- 1 1 ;t), Vn G N. (191) 

Consequently, for any n G N, we have: 

»%, = »X(0) = Z„({°°C*, °^*}jjli; 0) = °%, (192) 



and 



—iHot ooy iHot 



e 



Z n e lHat . (193) 



This completes our proof. 
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